Landau theory of glassy dynamics 
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An exact solution of a Landau model of an order-disorder transition with activated critical dy- 
namics is presented. The model describes a funnel-shaped topography of the order parameter space 
in which the number of energy lowering trajectories rapidly diminishes as the ordered ground-state 
is approached. This leads to an asymmetry in the effective transition rates which results in a non- 
exponential relaxation of the order-parameter fluctuations and a Vogel-Fulcher-Tammann divergence 
of the relaxation times, typical of a glass transition. We argue that the Landau model provides a 
general framework for studying glassy dynamics in a variety of systems. 



Introduction Glassy dynamics occur in a large variety 
of systems, such as supercooled liquids, foams and granu- 
lar matter. They are characterized by an extremely rapid 
increase of relaxation times and by a non-exponential de- 
cay of time-dependent correlation functions l'| . The rapid 
increase in time scales is typically fit by an exponential, 
Vogel-Fulcher-Tammann (VFT) dependence on a control 
parameter such as temperature or density. In this letter, 
we propose a Landau model of order-parameter dynam- 
ics in the vicinity of a critical point, which exhibits these 
features. 

The canonical theoretical framework for dynamics near 
a critical point is the time-dependent Ginzburg-Landau 
equation It successfully describes the phenomenon 
of critical slowing down, whereby the relaxation time 
for order-parameter fluctuations scales as a power of the 
correlation length, which in turn diverges as the critical 
point is approached. A natural question to ask is whether 
this framework can be adapted to describe glassy dynam- 
ics. Here we show analytically that a master equation, 
based on a Landau free energy, where the energy and 
entropy play asymmetric roles in determining the tran- 
sition rates, naturally leads to a VFT divergence and a 
broad distribution of relaxation time scales near a critical 
point. 

The microscopic basis for the asymmetry in the tran- 
sition rates can be found in the topology of the space of 
trajectories. Namely, we consider the situation where the 
number of energy lowering trajectories diminishes as the 
ordered ground state is approached. A similar funnel- 
like topology of the trajectory space for a simple kinetic 
model of protein folding has been identified recently as a 
possible explanation of fast folding'^. Here we show that 
such a topology of the trajectory space, when a system 
is driven to a critical point, can result in glassy criti- 
cal dynamics. Furthermore, we hypothesize that diverse 
systems such as foams, granular matter and supercooled 
liquids owe their glassy features to such a mechanism, 
and propose computational tests of this idea. 



Model We start with a mean-field model of a continu- 
ous phase transition for a scalar order parameter p which 
takes values on the integers. We describe dynamics of p 
as a random walk on the order parameter space which is 
a one-dimensional lattice. The Landau free energy in the 
disordered phase (/i < /j,*), is assumed to be quadratic 
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p is the dimensionless control parameter, and p* is its 
critical value. The energy and the entropy are also as- 
sumed to be quadratic functions of the order parameter: 
l3E{p) = (3E° - pp^/2, and S{p) ^ S° - p*p^/2 A\. The 
free energy f3F — j3E — S leads to a mean-field phase 
transition as p* — p 0. Below we shall be interested 
in the order parameter dynamics as the phase transition 
point is approached from the disordered phase. 

In the Landau approach to critical dynamics, P{p,t), 
the probability for the order parameter to take the value 
p at time t, is the solution to the master equation 



dP{p) 
dt 



= -{Wp^p-i + Wp^p+,)P{p) 



oP{p+^) -(2) 



The transition rates, Wp—,pi, depend on the free energy, 
the canonical choice being, Wp^p' = exp[— (/3_F(p') — 
/3F{p))/2)]. These rates satisfy detailed balance and 
yield the equilibrium Gibbs distribution, Pcq{p) oc 
exp[—PF{p)], as the steady state solution to Eq.j^Jl. Re- 
laxation processes described by these rates implicitly as- 
sume that neither energy nor entropy barriers impede 
the fluctuations of the order parameter around its equi- 
librium value, — 0. 

The canonical scenario described above can be altered 
dramatically if the microscopic dynamics are constrained 
in some way. An example is provided by a system of hard 
spheres at densities approaching random close packing. 
At these large densities only correlated motions of many 
spheres are available to the system in order to relax a 
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Scaling solution at the critical point For e = 0, the 



Laplace transform P{p, s) 
exactly by, 



P{p,t)e ^^dt is given 
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.M*|p|/2. 
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The long-time limit t oo corresponds to the s —^ 
limit in Laplace space and Eq.lO suggests the scaling 
limit: s — > and \p\ oo while keeping ^/se^ ''''/^ fixed. 



Using the scaling variable z — 



p|/2 



Vt 



, the scaling form 



of the solution to Eq.© can be written as 



FIG. 1: The Landau free energy, /3F(p) is a 
parabola whose width increases as the transition is 
approached (e — > 0). The order parameter p per- 
forms a random walk with rates dictated by the 
energy and the entropy, which are functions of p. 

density perturbation 5] . We consider the case where such 
constraints in the microscopic dynamics can be encoded 
as an asymmetry in the way energy and entropy enter 
the transition rates Wp^pi . In particular, we assume 
that energy lowering transitions can proceed only from a 
small subset of microstates at given p. If the system finds 
itself in one of these microstates the transition occurs at 
a fixed rate independent of the value of p. Thus the rate 
of energy lowering transitions is a function of the entropy 
change only. Taking the requirement of detailed balance 
also into account leads to 

I^p^p' = e-[^(^)-^^(''')l ■,Eip')<E{p) 

Wp^p. = e-['3^^(''')-^^('')l ■E{p')>E{p) . (3) 

Using the quadratic expressions for the energy and en- 
tropy in these transition rates leads to 

dP{p) 



at 



+ e-''''Pp{p-l) + e-''^P+'^^P{p+l) (4) 



for the master equation for P{p, t) for p > {). Since the 
rates are symmetric, the equation for p < has a similar 
form. Close to the critical point [p* — p <^ \) typical 
values of p are much greater than one, and we can take 
the continuum-p limit of this equation to obtain]^: 

dP 

'at 



— [epe^t^'Mp^on.-t-'lpl^^^ 



(5) 



dp ' dp ^ 

with e = p* ~ p the distance to the transition point, 
and D = 1/2. We note that a direct consequence of the 
asymmetric transition rates is the exponentially decaying 
factor multiplying both the drift and the diffusion terms 
in Eq. © . The time-invariant probability distribution ob- 
tained from this equation is the normalized, equilibrium 
distribution: Peq{p) - 



2\ogDt 



F{z) 



(7) 



The scaling function F{z) = ^e^^ is determined 

by substituting the scaling form, Eq.jTJ), in Eq.® with 
e = 0, and matching the leading order terms for large 
t. This solution immediately gives us the mean square 
fiuctuation of p for large i. 
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\og'{2p*'Dt). (8) 



The fluctuations grow logarithmically with time, in con- 
trast to the standard \/t scaling expected at a normal 
critical point. 

Scaling solution for e ^ Away from the critical 
point, the Fokker-Planck equation, Eq.©, has a nor- 
malizable stationary solution. The equation obtained by 
substituting P(p,i) = e-'-P^ ^'^^ip{p,t) in Eq.© allows a 
scaling solution for ip(p,t) in the limit < — > oo, p — > cx) 
but keeping the combination e'^ /t fixed. Repeating the 
same steps as described above for the e = case we com- 
pute 



PiP,t) 



exp {-ep^/AD 



2j^dp exp(-ep2/4£' 



(9) 



In the limit t — > oo, this solution recovers the equi- 
librium distribution. From Eq.lO, it is clear that at a 
given large there are two length scales in p-space: the 
equilibrium length scale, -^Z {p'^)eq = ■\/2_D/e, and a time- 
dependent length scale, pc — [1/ p*)\og{p*^Dt). For 
p <C Pc PiPjt) reduces to the equilibrium Gaussian dis- 
tribution while for p ^ pc the distribution of p is super- 
exponential. At long times pc becomes much larger than 
the equilibrium p and the system "knows" that it is in 
equilibrium. This defines a crossover time tc e^ v^W^ 
obtained by equating pc with \/ {p^)eq- Below we also 
compute the relaxation time from the the mean-square 
fluctuation of the order parameter, and show that 

it diverges as exp[l/e], which is precisely the VFT law. 
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FIG. 2: The probability distribution, P{p, t) for three 
different times, with e — 0.1, D — 1., and /i* = 
1. For these parameters, pc — 0.69, 2.3, 4.6 for 
t = 2, 10, 100, respectively, and = 88. Note 
that a,t t — 2 the distribution is essentiaUy super- 
exponential and for t — 100 it is a Gaussian. 



An interesting conclusion that can be drawn about the 
Landau model is that its critical dynamics cannot be 
characterized by a single diverging time scale. 

We have seen earlier that in the e = case, ~ 
log^(i) for large t. For e 7^ this quantity will relax to its 
equilibrium value {p'^)eq = 2_D/e. From Eq.®, we find 
in the scaling limit, 



Jo°°e.pi-epy4D-^^)dp ' 
Thus one can write, 

d\og[Z{A)] 



(10) 



dA 



(11) 



where A = e/AD and the generating function Z{A) is 
given by, 



Z[A) = / 
Jo 



dp exp ( —Ap 



oM P 



2p*^Dt 



From Ea. (|ll|l we find that for large t and small A, 
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(12) 



(13) 



Thus, at late times, the mean square fiuctuation relaxes 
to its equilibrium value in a power law ('^ 1/t) fashion. 
A standard way to estimate the relaxation time r is to 
define it as the time needed for the difference (p^)eq — 



(p^(t)) to decay to a given value of order one. From 
Ea. (|13() we immediately obtain 

D 



T OC 



(14) 



which diverges in the VFT manner as e — )■ 0. 

Activated dynamics at a critical point were previously 
shown to occur in systems with quenched disorder and 
it had been suggested that a similar situation rnight ex- 
ist in non-disordered systems with frustrationjSj. The 
emergence of a glass state in the latter case could then 
be associated with a critical point at which relaxation 
times would acquire astronomical values long before any 
discernible spatial correlations indicative of an ordered 
state could be established. The dynamical Landau model 
presented here provides an explicit realization of this sce- 
nario. 

The scaling solution to our model demonstrates that a 
seemingly innocuous change in the transition rates leads 
to an exponential divergence of timescales at a critical 
point, while all its static properties are described by 
mean-field Landau theory. The dynamics of the Lan- 
dau model are characterized by a multitude of timescales 
and by non-exponential relaxations even away from the 
critical point. Namely, the mean-square fluctuation of p 
can be written in terms of the density of states U{E) as: 
{p'^{t)) = / dEQ.{E)e-^K The power law decay for e > 0, 
Ea. (|13|) . implies that U,{E) is proportional to 5(£') — const 
with the ^-function guaranteeing the equilibrium value. 
At the critical point, the logarithmic decay of the fluctu- 
ations implies Vl{E) oc \t?{E)/E. The glass transition in 
our model is, therefore, not characterized by a vanishing 
gap in the energy spectrum (of Eq. ) but instead by an 
accumulation of states near E = Q. This then translates 
to a distribution of timescales with increasing weight in 
the tails of the distribution corresponding to long times. 

In the Landau model with a quadratic entropy func- 
tion the entropy goes to zero at pl^ax = 25'o//i* and our 
analysis is, therefore, strictly valid for {p'^)eq < pinax O'" 
e > p* /2Sq. For smaller values of e the relaxation time- 
scale divergence is cut off by the system size. 

Microscopic Models The transition rates given by 
Eq.Q, which formed the basis of our dynamical equa- 
tion and led to the VFT divergence of timescales, are 
explicitly realized in the interacting three-color modelj^. 
The configurations of the three color model on the hon- 
eycomb lattice can be mapped on to loop packings and 
the loops can be identified with equal height contours of 
a discrete height field which lives on the dual triangular 
latticejl^. Coloring configurations can be organized into 
topologically distinct sectors characterized by the num- 
ber of non-zero winding number loops, which correspond 
to different global tilts in the height representation. The 
introduction of a long-range interaction between one of 
the colors leads to a phase transition from a flat to a 



pnase 

completely tilted statep with the tilt playing the role of 
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the order parameter p. 

Microscopic dynamics for the three color model can 
be defined in terms of loop updates |^ with transitions 
between different loop configurations satisfying detailed 
balance. Simulations were performed to extract the tran- 
sition rates between different tilt states, Wp_p' and it was 
shown, numerically, that they take the form described 
by Eq.lPJ. The asymmetry in the rates owes its origin 
to the loop dynamics. A p-reducing (energy-increasing) 
transition involves deleting an existing non-zero winding 
number loop. The rate of this transition is determined 
only by the change in energy since all loop configurations 
at a given p can lead to a configuration with p' < p. 
On the other hand, a transition that increases p and 
consequently decreases the energy, involves introducing 
a new non-zero winding number loop. Since only a small 
subset of loop configurations at a given p can accommo- 
date a new non-zero winding number loop, this transition 
rate depends only on the change of entropy. We believe 
similar mechanisms are operative in other system where 
transitions between macro-states involve large scale rear- 
rangements which, due to the presence of constraints, can 
occur only through coordination of many microscopically 
rearranging regions. 

In this work, we identify a concrete physical mecha- 
nism by which dynamical barriers are established in the 
free energy landscape. The dynamics in the order pa- 
rameter space in our model is thus reminiscent of diffu- 
sion on hierarchical lattices^,^^ ^^'^ trap models^lSj]. 
Recently, attempts have been made to identify inher- 
ent structures in supercooled liq uids as providing a mi- 
croscopic basis for traps^^^ llSj . If we identify inher- 
ent structures with the order parameter p in our model 
then the predicted distribution of trapping times, Tp = 
l/{Wp^p-i + Wp^p+i), is log-normal with a width that 
increases as 1/e . This is consistent with the behavior 
observed in supercooled liquid simulations l^]. 

Studies of kinetically constrained models [IS] have iden- 
tified dynamical heterogeneities [tR resulting from the 
kinetic constraints and an entropy crisis in trajectory 
space which leads to rapidly increasing timescaleS (18i.il9l] . 
In our model, the asymmetry of the rates results in an 
entropy crisis and, combined with the existence of an 
equilibrium critical point, leads to an exponential diver- 
gence of relaxation time-scales for a finite value of the 
control parameter. 

The coarse-grained dynamical model presented here 
suggests a new way of analyzing atomistic simulations of 
glassy materials. From molecular dynamics simulations 



one can measure the transition rates between macro- 
states characterized by different values of a slow variable 
which plays the role of an order parameter, for exam- 
ple, the metabasin energies in supercooled liquids 
This approach avoids making specific assumptions about 
the origin of the transition rates. Dynamics in order pa- 
rameter space can be constructed explicitly and analyzed 
for the types of asymmetries and related dynamical bar- 
riers discussed here within a simple Landau model. Work 
along these lines is currently in progress. 
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